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I.  INTRODUCTION 


The  dependence  of  the  threshold  voltage  and  radiation  response  of 
n-channel  AlGaAs/GaAs  modulation  doped  field-effect  transistors  (MODFETs) 
on  acceptor  doping  density  has  been  analyzed  previously  (Refs.  1  and  2). 
These  analyses  have  been  extended  to  describe  the  dependence  of  MODFET  I-V 
characteristics  on  acceptor  doping  density.  A  triangular-well,  one- 
subband,  depletion  layer  model  has  been  developed  that  applies  over  the 
range  of  I-V  characteristics  from  subthreshold  to  saturation,  some  nine 
orders  of  magnitude  in  drain-source  current. 

1 1 

For  typical  unintentional  acceptor  doping  densities  of  10  to 
10^^  cm“^,  characteristic  in  molecular  beam  epitaxy  (MBE)  grown  structures, 
we  show  that  the  experimentally  derived  threshold  voltage  differs  from  the 
strong  inversion  model  threshold  voltage  (Ref.  1)  by  0.25  V  at  jceptor 
densities  of  10^^  cm~^.  At  acceptor  densities  of  10^^  cm~^,  the  difference 
between  the  strong  inversion  model  and  the  experimental  extrapolation  for 
the  threshold  voltage  is  about  0.12  V. 

Inclusion  of  the  acceptor  doping  density  is  shown  to  account  for  the 
discrepancy  between  the  AlGaAs  layer  and  the  device  capacitance  per  unit 
area  described  in  the  literature  (Ref.  3). 
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II.  MODEL 


The  band  structure  of  a  typical  AlGaAs(n)/GaAs  hetero junction  with  a 
Schottky  barrier,  4)^,  at  the  gate  and  a  spacer  layer  at  the  interface  under 
bias,  Vg,  is  shown  in  Fig.  1. 

Under  the  restrictions  imposed  by  the  assumptions  cited  in  the  caption 
of  Fig.  1,  Poisson's  equation  may  be  integrated  across  the  structure  to 
yield  the  applied  gate  voltage  as  a  function  of  device  geometry,  doping 
densities,  and  channel  charge,  n^: 

Vg^V^.fCn^)  (1) 

where  is  the  difference  between  the  Schottky  barrier  height  and  the  sum 
of  the  AlGaAs/GaAs  band  offset  and  potential  drop  across  the  doped  AlGaAs 
layer  resulting  from  the  ionized  donors.  The  function  f(ng)  may  be  written 
as 


f(n  )  =  (q/e)(d  +  a)(N  W  +  n^)  +  C  (N  W  +  +  { kT/q )ln[ exp( n  n  )  -  1] 

(2) 

where  is  a  function  of  the  Planck  constant,  the  carrier  effective  mass,  the 

elemental  charge,  and  the  permittivity  of  AlGaAs  and  GaAs,  all  assumed  to  be 

equal.  is  equal  to  ~  1.7  x  10  ^  V-cm^'^^.  Similarly,  the  charge  density  n^, 

is  a  function  of  physical  constants  and  the  effective  mass  of  the  carriers  and 

1  1  P 

is  equal  to  ~  8.4  x  10  cm"  .  In  the  next  section,  we  will  exploit  the 
mathematical  properties  of  the  function  f(ng)  to  derive  the  electrical  proper¬ 
ties  of  these  devices. 
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Fig.  1.  Band  Diagram  of  a  Typical  AlGaAs(n)/GaAs  MODFET  with  Schottky  Gate, 

Under  Bias.  In  the  depletion  layer  approximation,  the  donors  and 

acceptors  are  assumed  to  be  completely  ionized  in  the  doped  AlGaAs 

layer,  d,  the  spacer  layer,  a,  and  the  depletion  layer,  W.  The 

doping  densities,  and  N^,  are  assumed  to  be  constant.  A  delta- 

function  channel  charge  distribution  at  the  average  channel  width  is 

assumed.  Band  bending  from  the  interface  at  (d  +  a)  to  the  edge  of 

the  depletion  region  (W  +  d  +  a)  is  the  difference  of  position  of 

conduction  band  relative  to  the  Fermi  level,  E  /2  +  41.  ,,  ,  and  the 

g  bulk 

Fermi  level  relative  to  the  bottom  of  the  two-dimensional  channel. 
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III.  ELECTRICAL  PROPEHTIES 


A.  THRESHOLD  VOLTAGE 

At  threshold,  we  require  that  the  channel  density  be  equal  to  the 

acceptor  density,  N^,  times  the  average  channel  width,  ,  which  may  be 

calculated  in  the  triangular-well  approximation  using  variational  wave 

functions  (Ref.  M).  This  definition  for  threshold  is  consistent  with  the 

strong  inversion  definition  of  threshold  in  metal  oxide  semiconductor  field 

effect  transistors  (MOSFETs).  Our  definition  is  the  two-dimensional 

equivalent.  Evaluating  Eq .  (2)  at  threshold  and  substituting  into  Eq.  (1) 

yields  the  threshold  voltage.  The  results  of  this  calculation  indicate  that 

the  threshold  voltage  is  very  sensitive  to  the  acceptor  doping  density  above 
1 4  T 

-  10  cm  The  details  of  this  calculation  have  been  presented  elsewhere 
(Ref.  1). 

B.  SUBTHRESHOLD  I-V  CHARACTERISTICS 

In  the  subthreshold  region,  where  n^  <  over  the  whole  channel, 

f(ng)  may  be  approximated  as  follows: 

f(ng)  =  (kT/q)lnlg(ng)]  (3) 

Solving  Eq.  (3)  for  gin^),  expanding  to  first  order  about  the  threshold 
charge,  and  using  Eq.  (2)  to  determine  the  Taylor  series  expansion  coeffi¬ 
cients  yield  a  convenient  approximation  for  f(ng)  when  substituted  back  into 

Eq .  (3).  Shown  in  Fig.  2  is  the  function  f(n  ),  which  is  equivalent  to  (V_  - 

^  6 

Vq)  vs  the  lug  of  the  channel  charge  for  two  values  of  the  acceptor  doping 
density.  The  long  dashes  are  the  results  of  the  approximation  just  described. 

Using  the  subthreshold  approximation  for  f(ng)  in  Eq.  (1),  and  solving 
for  n  in  terms  of  \/„,  we  may  calculate  the  subthreshold  I-V  characteristics 
using  a  charge  control  model  (Ref.  5).  The  results  of  this  calculation  yield 
the  MODFET  equivalent  of  the  MOSFET  charge  sheet  subthreshold  characteristics 
(Ref.  6).  Details  of  these  results  will  be  discussed  at  the  end  of  this 
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Vo)  [Volts] 


Channel  Charge  Density  [cm 

Fig.  2.  Vg  -  vs  the  Log  of  ♦'he  Channel  Charge.  The  solid  lines  are  the 

results  given  by  Eq.  (2).  The  long  dashes  are  the  results  of  the 
Taylor  expansion  of  f(ng)  in  subthreshold,  as  described  in  the 
text.  For  an  acceptor  doping  density  of  10^^  cm"^,  the  threshold 
charge  density  is  less  than  10°  cm~^.  Therefore,  the  subthreshold 
expansion  is  useful  well  above  the  threshold  charge  density.  At  an 
acceptor  doping  density  of  10^^  cm‘^,  the  threshold  density  is  near 
10^^  cm”^,  and  this  expansion  is  only  useful  for  channel  charge 
densities  below  threshold.  The  short  dashes  are  the  results  of  a 
Taylor  expansion  of  f(ng)  about  n^,,  as  described  in  the  text. 
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see^ iOn  in  the  context  of  the  complete  desct-iption  of  I-V  chafacterist ics  from 
^aubthreaho  id  to  saturation. 

C .  SATUHATION  I-V  CHARACTERISTICS 

We  del'ine  tiie  saturation  region  such  that  n^  •  n^  over  the  whole  channel. 

In  this  regiorj,  Eq .  (c)  may  be  expanded  in  a  Taylor  series  in  n^  about  n^. 

The  results  of  this  approximation  to  first  order  are  sriown  as  short  dashes  in 

Fig.  2.  Above  le,  (8.4  «  10^^  cm""^ ) ,  the  expar.sion  is  quite  good.  Much  below 
-2 

n^,  (  d  .  Ill  cm  ),  the  first  order  expansion  depar'ts  fr-om  the  exact  result 

and  apprua.mris  a  constant  at  low  channel  densities  for  all  acceptor  densities. 

As  bei^.au;,  the  approximation  for  f(n^)  may  be  substituted  into  Eq.  (1),  which 

is  invcrteu  to  yield  r.  ,  as  a  function  of  V, .  The  result  is 

^  g 

n  :  n  +  K  ^[V  -  V  -  f ( n  )  j  '  ( kT  q )  (4) 

sc  g  o  c 

wher'tj  K  is  a  constarit  that  depends  on  the  device  geometry,  doping  densities, 
depletion  width,  and  physical  constants.  This  form  for  Is  different  than 
the  form  previously  assumed  (Hef.  5).  The  previous  form  ignores  the  contri¬ 
bution  from  ri  and  f(n„)  and  implicitly  assumes  that  is  (kT,'q)C.,„  .  . 

^  C  M  2.uiiAS 

The  charge  control  model  may  then  be  used  to  calculate  saturation  1-V 
characteristics.  We  will  defer  the  detailed  discussion  of  these  results  to 
the  end  of  this  section. 

D.  DEVICE  CAFACITANCE 

The  device  capacitance,  when  the  charge  density  in  the  channel  is  greater 

than  n  ,  may  be  determined  by  differentiating  Eq .  (4)  with  respect  to  V  .  The 
^  6 

resulting  capacitance  per  unit  area,  may  be  cast  in  the  following  form 

al  “ci 

(Ref.  3): 


C  -c/(d-t-a  +  Ad) 
area 


where  Ad  is  given  by 


Ad 


(2t/3q)C  (N  W  +  n  ) ^ +  1 .58( kTc/q^ ) /n 
o  a  c  c 


(5) 


(6) 


Whei,  N  W  ri  , .  Ad  •  wl-  A.  Fot  itii'ge  valLitri  of  th*;  acceptor  dup;!.c  Jri.oji, 
ly'’*  -  s  *’ 

(  -  it’  cn.  I,  Ad  1:3  r'.rd..ctd  tu  -  75  A.  The  coriritant  K,  in  the  pi'c-viocc  jct- 
cection,  1  i  related  t'^  t n*:.  capacitaricc  p’e’i  jnit  area  a;:>  toiirii^ai 


K 


(kT. 


)C. 


dlfOa 


E .  E X h£h i MEN! AL  THhESHOLh  VOLTAGE 

MijIiFEl  trc'esholr!  vfdtaget;  are  determined  exper  inict;ta.  ly  oy  e>;  t  rapoi  a  t  i  ng 
tne  aataration  curr'ent,  or  squar-e  root  of  the  saturation  current,  vs  gate 
voltage  to  zer'O.  The  gate  voltage  intercept  is  the  exper imer.taJ  ly  deter  mined 
threshold  voltage.  This  is  mathematically  equivalent  to  solving  Eq.  (4)  for 
the  gate  voltage  when  ti^  is  equal  to  zero.  This  calculation  yields  an  experi- 
mental  threshold  voltage  that  differs  from  the  strong  niversion  definition  of 
the  tnreshiild  voltage.  The  difference  between  the  experimental  vaiuc  and  the 
theoretical  value  is  giver  by 


wher'e  f(h(.j|') 
density.  In 
the  acceptor 


•tn  =  -  ^"c  "-'area 

is  the  vaoue  of  Eq.  (2)  when  evaluated  at  the  thresnoid  charge 
Fig.  3,  the  threshold  voltage  difference  is  plotted  vs  the  log  of' 
detisi  ty . 
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Threshold  Voltage  Difference  (Vex  ''^th) 

[Volts] 


Fig.  3.  Difference  of  the  Experimental  and  the  Strong  Inversion  Threshold 

Voltage  vs  the  Log  of  the  Acceptor  Density.  Each  term  in  Eq.  (8)  is 
plotted  separately.  At  low  acceptor  densities  (<  10^^  cm“^),  the 
difference  may  be  as  much  as  0.25  V.  This  difference  decreases  as 
the  acceptor  density  increases.  Because  this  difference  depends  on 
the  acceptor  density,  a  comparison  of  experimental  threshold  voltages 
may  not  be  appropriate  if  acceptor  densities  differ  significantly. 
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IV.  I-V  CHARACTERISTICS  FROM  SUBTHRESHOLD  TO  SATURATION 


In  the  charge  control  model,  1-V  characteristics  are  determined  by 
substituting  V_  -  V-(x)  for  V„  in  Eq.  (1),  inverting  the  result  to  find  n„  as 

6^6  i> 

a  function  of  Vj,{x),  substituting  the  result  into  the  relationship  for  the 
current  at  position  x  in  the  channel,  and  integrating  over  the  channel  length 
(Ref.  5).  In  the  subthreshold  and  saturation  regions,  as  defined  in 
subsections  III.B  and  III.C,  the  approximations  for  fCng)  permit  a 
straightforward  inversion  of  Eq.  (1)  for  this  purpose.  As  seen  in  Fig.  2, 
there  is  a  region  over  which  neither  the  subthreshold  nor  the  saturation 
expansion  applies.  Because  the  integration  for  the  current  is  performed  over 
the  voltage  in  the  channel  as  well  as  over  the  channel  length,  an 
approximation  must  be  determined  that  is  continuous  in  the  voltage,  may  be 
inverted  to  explicitly  determine  the  channel  charge  in  terms  of  the  channel 
voltage,  and  reasonably  approximates  f(ng)  over  the  region. 

To  satisfy  these  requirements,  we  have  derived  a  piecewise  approximation 
for  ('(Og)  over  this  region.  From  the  threshold  charge  to  the  channel  charge, 
n^,  at  which  the  exact  ("(rig)  is  halfway  between  the  subthreshold  approximation 
at  threshold  and  the  saturation  approximation  at  n^,  we  have  approximated 
f(n„)  as  the  log  of  a  linear  function  of  n_.  The  two  expansion  coefficients 
are  chosen  so  that  this  approximation  for  f(n„)  is  exact  at  n„  and  connects 
with  the  subthreshold  approximation  at  the  threshold  channel  density.  From  n^ 
to  n^,  we  have  assumed  that  f(ng)  is  linear  in  ng.  The  two  expansion  coeffi¬ 
cients  are  determined  so  that  the  approximation  in  this  region  is  exact  at  the 
extremes,  n^  and  n^,. 

The  charge  control  model  may  now  be  used  to  calculate  the  I-V  character¬ 
istics.  A  complication  arises  in  the  application  of  the  charge  control  model 
because,  for  various  values  of  the  applied  gate  and  drain-source  voltages, 
different  regions  of  the  channel  may  have  charge  densities  that  must  be  calcu¬ 
lated  by  different  approximations  to  f(ng).  Therefore,  the  current  equation 
must  be  integrated  placewise,  although  never  over  more  than  four  regions,  as 
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our  piecewise  approximation  for  fCn^)  requires.  Over  a  range  of  gate  (or 
drain-source)  voltages,  10  possible  combinations  of  regions  might  occur. 
Shown  in  Fig.  4  is  the  drain-source  current  vs  drain-source  voltage,  for  the 
parameters  given,  for  various  gate  voltages  using  the  approximations  Just 
described. 

In  Fig.  5,  we  show  the  drain-source  current  vs  applied  gate  voltage  for 
three  values  of  the  acceptor  density  id  two  values  of  the  drain-source 
voltage. 

Intrinsic  transconductances  may  be  calculated  by  differentiating  the 
drain-source  current  relationships  with  respect  to  the  gate  voltage.  Having 
done  this  calculation,  we  obtain  transconductances  of  ~  300  mS/tnm  for  struc¬ 
tures  with  acceptor  doping  of  10^^  cm“^,  evaluated  at  zero  gate  voltage  and 
2.5  V  drain-source  voltage. 
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Drain— Source  Current  [mA] 


Fig. 


4.  Drain-Source  Current  vs  Drain-Source  Voltage.  The  current  saturates 
at  large  drain-source  voltage  without  velocity  saturation  or  cutoff 
of  the  model  invoked.  When  the  whole  channel  is  forced  into  the 
sub threshold  charge  region,  the  current  depends  on  a  constant  plus  an 
exponential  in  -V^g  which  is  negligible  at  large  drain-source 
voltages. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories; 

Acrophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  prt^llant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-fleld-of-view  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistiy. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometiy,  micro- 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  if  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
eflects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


